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the nondegenerate case. In this work, we further study the pair and associated production of the 
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enhancement from cascade decays of the charged scalars H^^, We consider three important 
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and 200 GeV. The neutral scalars in the negative scenario are thus accessible at LHC run II. 
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I. INTRODUCTION 

In a previous paper |[T1, we presented a comprehensive analysis on the LHC signatures of the type II 
seesaw model of neutrino masses in the nondegenerate case of the triplet scalars. In this companion paper, 
another important signature—the pair and associated production of the neutral scalars-is explored in great 
detail. This is correlated to the pair production of the standard model (SM) Higgs boson, h, which has 
attracted lots of theoretical and experimental interest [3l since its discovery lH |5l , because the pair 
production can be used to gain information on the electroweak symmetry breaking sector Q. Since any 
new ingredients in the scalar sector can potentially alter the production and decay properties of the Higgs 
boson, a thorough examination of the properties offers a diagnostic tool to physics effects beyond the SM. 
The Higgs boson pair production has been well studied for collider phenomenology in the framework of the 
SM and beyond Il6l422l. and extensively studied in various new physics models ll23 - l47i . as well as in the 
effective field theory approach of anomalous couplings lldStiSTTl and effective operators |[52] - [56]l . 

The pair production of the SM Higgs boson proceeds dominantly through the gluon fusion process 10|8l, 
and has a cross section at the 14 TeV LHC (LHC 14) of about 18 fb at leading order f^. [^It can be utilized 
to measure the Higgs trilinear coupling. A series of studies have surveyed its observability in the 6677 , 
hhW~^W~, bbbb, and WW*WW* signal channels ifTTMTSl 1211 l22l 1^l57l . For the theoretical 
and experimental status of the Higgs trilinear coupling and pair production at the LHC, see Refs. 157115^ . 
In summary, at the 14 TeV LHC with an integrated luminosity of 3000 fb“^ (LHC14@3000), the trilinear 
coupling could be measured at an accuracy of ~ 40% IT^ . and thus leaves potential space for new physics. 

As we pointed out in Ref. lT|, in the negative scenario of the type II seesaw model where the doubly 
charged scalars are the heaviest and the neutral ones H^/A^ the lightest, i.e., Mh±± > Mh± > 
MjjojAQ, the associated production gives the same signals as the SM Higgs pair production while 

enjoying a larger cross section. The leading production channel is the Drell-Yan process pp —)• Z* —)• 
H^A^, with a typical cross section 20-500 fb in the mass region 130-300 GeV. Additionally, there exists a 
sizable enhancement from the cascade decays of the heavier charged scalars, which also gives some indirect 
evidence for these particles. The purpose of this paper is to examine the importance of the H^A^ production 
with an emphasis on the contribution from cascade decays and to explore their observability. 

The paper is organized as follows. In Sec.|^ we summarize the relevant part of the type II seesaw and 
explore the decay properties of , A^ in the negative scenario. Sections[^and[^contain our systematical 
analysis of the impact of cascade decays on the H^/A^ production in the three signal channels, 6677 , 
bbT^T~, and bbi^£~^T- We discuss the observability of the signals and estimate the required integrated 


' This number is modified to 33 fb at next-to-leading order Q and to 40 fb at next-to-next-to-leading order [161. 
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luminosity for a certain mass reach and significance. Discussions and conclusions are presented in Sec. 
In most cases, we will follow the notations and conventions in Ref. 


II. DECAY PROPERTIES OE NEUTRAL SCALARS IN THE NEGATIVE SCENARIO 


The type II seesaw and its various experimental constraints have been reviewed in our previous work 
[Tj. Here we recall the most relevant content that is necessary for our study of the decay properties of the 
scalars in this section and of their detection at the LHC in later sections. 

The type II seesaw model introduces an extra scalar triplet A of hypercharge two |[59l on top of the SM 
Higgs doublet of hypercharge unity. Writing A in matrix form, the most general scalar potential is 

L(^>, A) = + M^Tr{A^A) + + A 2 (Tr(AtA))^ + A3Tr(Al'A)2 

+A 4 (T>t«>)Tr(A^A) + A 5 T>tAAtT>+ + H.c.) . (1) 


As in the SM, m? < 0 is assumed to trigger spontaneous symmetry breaking, while > 0 sets the mass 
scale of the new scalars. The vacuum expectation value (vev) r; of <h then induces via the /x term a vev va. 
for A. The components of equal charge (and also of identical CP in the case of neutral components) in 
A and then mix into physical scalars A^-, , h and would-be Goldstone bosons with the 

mixing angles specified by (see, for insfance. Refs. ll^lMI ) 


tan 0 + = 


V 


tana = 


2va 

V 


tan 200 


2 va n^(A4 + As) — 2M^ 

V 2v‘^\i-M\-v\{\2 +X'iY 


( 2 ) 


where an auxiliary parameter is introduced for convenience. 


Mi = 


'J^va 


(3) 


To a good approximation, the SM-like Higgs boson h has the mass ~ 'j2X[v, the new neutral scalars 
have an equal mass M^o k, M^q k, Ma, and the new scalars of various charges are equidistant in 
squared masses: 


M 




-M 




Mfj± - M 


HO/AO 




(4) 


There are thus two scenarios of spectra, positive or negative, according to the sign of A 5 . For convenience, 
we define AM = M^± — M^o 

In the rest of this section, we discuss the decay properties of the new scalars in the negative scenario 
with an emphasis on and A^. The explicit expressions for the relevant decay widths can be found in 
Refs. Il62l - l64l . It has been shown that /A^ decays dominantly into neutrinos for va < 10“^ GeV ll65l . 
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FIG. 1. Branching ratios of and versus Ma at some benchmark points of AM and va- {AM,va) 
(5,0.01), (10,0.01), (5, 0.001) GeV, from the upper to the lower panels. 





















5 


resulting in totally invisible final states. We will restrict ourselves to va S> 10“^ GeV in this work, 
where H^/A^ dominantly decays into visible particles. Before we detail their decay properties, we give 
a brief account of the cascade decays of the charged scalars. The branching ratios of the cascade decays 
are controlled by the three parameters, va, AM, and Ma- The cascade decays dominate in the moderate 
region of va and for AM not too small, where a minimum value of AM ~ 2 GeV appears around 
VA ~ 10“"^ GeV im[6211^1^. In Fig. the branching ratios of and are shown as a function 
of Ma at some benchmark points of va and AM. Basically speaking, in the mass region Ma = 130- 
300 GeV, the cascade decays are dominant for a relatively large mass splitting AM (as shown in the 
middle panel of Fig. [TJ or a relatively small va (in the lower panel). 

A. decays 

At tree level, can decay to // (/ = q, 1), vv, W^W~, ZZ, and hh. It can also decay to gg, 77 , 
and Z 7 through radiative effects. Similarly, //, uv, Zh at tree level, and it has the same decay 

modes as at the loop level. Since we have chosen va ^ 10“^ GeV, the neutrino mode can be safely 
neglected for both and A^. Previous work usually concentrated on the decoupling region where the 
neutral scalars H^/A^ are much heaver than the light CP-even Higgs h and the scalar self-couplings A* are 
taken to be zero for simplicity ll65l . In this case, the mixing angle Oq ^ a, and the H^W^W~ coupling 
[being proportional to sin(a — 0o)] tends to vanish. As a consequence, the VF-pair mode is absent and the 
dominant channels are —)• hh, ZZ for a heavy H^. In contrast, we take into account the effect of scalar 

self-interactions and focus on the nondecoupling regime, i.e., H^/A^ are not much heavier than h. 

For illustration, we choose the benchmark values va = 10“^ GeV, AM = 5 GeV; then, A 5 is deter¬ 
mined by Eq. Q upon specifying Ma- To investigate the effect of the scalar self-interactions, we note 
the following features in the decays of 1) The decay widths of —>• //, gg differ from those of h 

only by a factor of sin^ Oq, which leads to similar behavior for and h. 2) The only free parameter for 
the mixing between and h is A 4 , because [as shown in Eq. (j^J the impact of A 2,3 is suppressed by a 
small va and a relatively large mass difference between Ma and M^ while Ai is fixed by M^. 3) A 4 enters 
the and H^ZZ couplings and thus affects the decays —)• W^W~, ZZ. 4) The H^hh 

coupling simplifies for ua n such fhaf fhe only free paramefer in fhe decay —)• hh is again A 4 . Asa 
consequence of fhese feafures, we shall choose A 4 as a free parameter and vary if in fhe range [—1.0,1.0], 
and fix fhe couplings A 2 = A 3 = 0.1 which are involved in loop-induced decays. 

We firsl examine fhe branching ratios of —)• //. BR(P° —)• bb) and BR(iT° —)• ti) are plotted in 


^ As pointed out in Ref. 1621 . varying va in the range 10 ^-1 GeV would not change the branching ratios significantly. 
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FIG. 2. Branching ratios of —?► bb and tt as a function of Mfjo for various values of A 4 . 

Fig. 1^ for different mass regions of It is clear that the variation of BR{H^ —)■ bb) is more dramatic 

for A 4 > 0. The maximum of BR{H^ —)• bb) appears at A 4 ss 0.5. Obviously, BR(iF’^ —)■ bb) is a 
nonmonotonic function of A 4 , while BR{H^ —)• tt) monotonically increases with A 4 . As will be discussed 
later, this different behavior in the two mass regions is due mainly to a zero in the H^ZZ coupling. 

Now we study the bosonic decays —)• W^W~, ZZ, hh. In the left panel of Fig. we present 
the branching ratios of —)■ W~^W~, ZZ in the mass region 130-300 GeV. For most values of A4, 
BR(iT0 ^ W+W- ■) increases with Mj^o when Mjjo < 2Mw, and varying A4 for A4 > 0 changes it 
considerably. A4 has a strong impact on BR(Ff° —)> W~^W~) in the mass region 2Mz < M^o < 2Mh 
where the decay channel dominates overwhelmingly for A4 < 0 but becomes negligible for A4 approaching 
about 0.5. However, once the —)• hh channel is opened, —)• W~^W~ is suppressed significantly 

independent of A4. The decay ZZ cannot dominate when < 2Mw- In the mass region 2Mz < 

Mfjo < 2Mh, it is complementary with the W~^W~ channel, so their behavior is just opposite. More 
interestingly, there is a zero point for the H^ZZ coupling, which is proportional to {v sin 9q — dwA cos 60 ). 
According to Eq. Q, one obtains the corresponding Ma at the zero: 

MiiZZ) = ^2Ml-^{X4 + X5W. (5) 

Note that the above relation only holds for A 4 + A 5 < ss 0.5, since we are working in the scenario 

where Ma > Mh- The existence of the zero coupling explains the presence of the nodes in BR(iF‘^ ZZ) 
for A 4 < 0. 

^ The influence of A 4 for light fermions 6 , c, r, /j, and gluons is similar, so we only present —>■ bb) in Fig. 2 
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FIG. 3. Left; Branching ratios of 77° — W'^W , ZZ as a function of M^o in the mass region 130-300 GeV. Right; 
Branching ratios of 77° —>■ hh, ZZ as a function of M^o in the mass region 200-1000 GeV. 

In the right panel of Fig. BR(77° —)■ hh, ZZ) are shown in the mass region 200-1000 GeV. When 
Mfjo > 2Mh, the dependence on A 4 is simple: a larger A 4 corresponds to a smaller BR(77° —)• hh) and a 
larger BR(77° —)• ZZ). It is clear that A 4 has a more significant impact in the mass region 200 ~ 350 GeV, 
and varying A 4 could change BR(77° —)• ZZ) from 0 to 0.9. Once Mj^o exceeds 2Mt, the evolution of 
Br(77° —)• hh, ZZ) becomes smooth with the increase of Mj^o. There also exists a zero point for the H^hh 
coupling, which can be obtained as for the ZZ channel: 


Ml{hh) = ^2(A4 + A5)n2 -2M2 , 


( 6 ) 





























FIG. 4. Branching ratios of H'^ 77 , Z 7 as a function of M^jo for various sets of A 2,4 values, 

which is valid for A 4 + A 5 > 3M|/2t;^ ss 0.375. 

Finally, we investigate the loop-induced decays, —)■ 77 , Z 7 . In addition to the usual contributions 

from the top quark and W boson, the new charged scalars and also contribute to the decays. 
These new terms involve the and H couplings, which are proportional to 

: [( 2 A 2 + 2 A 3 — As) sin a cos 0o ~ ( 2 A 4 + A 5 ) cos a sin 0o], 

: (A 2 sin a cos 6*0 — A 4 cos a sin 9o). (7) 

One therefore has to consider the scalar self-couplings A 2 , 3 . For simplicity, we set A 2 = A 3 and vary them 
from —3.0 to 3.0. In Fig. we display BR(F7° —)> 77 ) and BR(F7° —)> Z 7 ) versus Mfjo for some typical 
sets of A 2,4 values. The evolution of both branching ratios crosses 3 orders of magnitude in this parameter 
region. The resulting enhancement compared with /i —;■ 77 in the SM looks significant: the maximal 
enhancement can be achieved at the level of 9% for the —)• 77 channel at Mfjo = 130 GeV, and of 

0.7% for the —)> Z 7 channel at M^o ~ 140 GeV. 

B. A° decays 

Similar to H^, the decay widths of —)• //, gg differ from those of /i by a factor of sin^ a with a 
being given in Eq. (|^. Moreover, the only vertex which involves A* is the A^Zh coupling proportional to 
(cos 9q sin a — 2 sin 9q cos a). As a consequence, one can only choose A 4 as a free parameter to illustrate 
the influence of scalar interactions. In this section, we also vary A 4 from —1.0 to 1.0 and take the same 
benchmark values for wa and AM as for the decays. 
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FIG. 5. Branching ratios of > bb, tt as a function of M^o for various values of A4. 


In the left panel of Fig. we present —)• hb) as a function of M^o. For a fixed value of A 4 , 

BR( 74 ‘^ —)• bb) decreases as M^o increases. The dependence of BR(A'^ —)• bb) on A 4 is simple: The larger 
A 4 is, the larger BR(^‘’ —)• bb) is. And BR(A° —)■ bb) can be dominant with A 4 = 1.0 as long as —)• Zh 
is not fully opened. The right panel of Fig. [^shows BR(A° —)• tt), which is very similar to BR{H^ —)• tt). 
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FIG. 6 . Branching ratios of > Zh as a function of M^o for various values of A4.. 

We then study the most important decay A^ —)• Zh. In Fig. we present BR(A‘^ —)• Zh) as a function 
As before, the influence of A 4 on the — >■ //, gg channels is similar to the bb mode. 
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FIG. 7. Branching ratios of > 77 and Z 7 as a function of M^o for various values of A 4 . 

of M^o in the low-mass region (130-300 GeV) and high-mass region (300-1000 GeV), respectively. The 
evolution of BR(74° —)• Zh) with M^o and A4 is just opposite to that of ^ 4 *^ —)• hb {tt) in the low- (high-) 
mass region. The variation of BR(A'^ —)• Zh) with A4 is dramatic below the Zh threshold. In particular, 
near the Zh threshold BR(74'^ —)• Zh) ~ 1.0 for A4 = —1.0, while BR(A‘^ —)• Zh) tends to vanish for 
A4 = 1.0, which corresponds to the zero point of the A^Zh coupling: 

Ml{Zh) = yj{\i + X^)v‘^ - Ml (8) 

with A4 -|- A5 > 0.5. BR(A'^ —)■ Zh) is totally dominant in the mass region between the Zh and 

tt thresholds, and becomes comparable to BR(74‘^ tt) when M^o > 2 Mt. 

At last, we study the one-loop-induced decays, —)■ 77 , Z 7 . These two channels can only be induced 
by the top quark in the loop since the A^W~^W~, A^H'^H~, and A^H^~^H couplings are absent in the 
GP-conserving case. In Fig. both BR(A° —)■ 77 ) and BR(A‘^ —)• Z 7 ) are displayed. For M^o below the 
Zh threshold, the variation in A4 of BR(A*^ —)• 77 ) increases as M^o increases. BR(A° —)• 77 ) could reach 
9 X 10“^ for M^o « 210 GeV and A 4 = 1.0, which is much smaller than the maximum of BR(F7° —)> 77 ). 
The variation in A 4 of BR(A° —)• Z 7 ) is slightly steeper, with a maximum of 1.2 x 10“^ at M^o 215 GeV 
and A4 = 1.0. 

In the above, we have discussed the decay channels of and A^ separately. We have shown that the 
scalar self-interactions have a large impact on their branching ratios. In Sec.|^ we will explore their LHC 
signatures. For this purpose, we choose the following benchmark values: 


UA = 0.001 GeV, AM = 5 GeV, A 2 = A 3 = 0.1, A 4 = 0.25. 


( 9 ) 
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FIG. 8. Branching ratios of /A^ as a function of M^o at the benchmark point in Eq. (j^. 


The reason that we set relatively small values of and AM is to obtain large cascade decays of charged 
scalars as well as a large enhancement of neutral scalar production. In Fig. we display all relevant 


branching ratios versus M^o for this benchmark model, which is to be simulated in Sec. IV for the LHC 
in the 6677 , bbT^T~, and bbW~^W~ signal channels. 


III. PRODUCTION OF NEUTRAL HIGGS FROM CASCADE DECAYS 

We pointed out in Ref. lUl the importance of the associated production in the nondegenerate case. 

To estimate the number of signal events, we simulated the signal channel at = 130 GeV. 

We found that, with a much higher production cross section than the SM Higgs pair {hh) production, a 
2.9cr excess in that signal channel is achievable for LHC 14@300. In the present work, we are interested 
in the observability of the associated production in the nondecoupling mass regime (130-200 GeV). 
In Fig. I^we first show the production cross sections for a pair of various scalars at LHC 14 versus Ma 
with a degenerate spectrum. As before, we incorporate the next-to-leading-order (NLO) QCD effects by 
multiplying a iC-factor of 1.3 in all qq production channels f7j|. The hh production through gluon-gluon 
fusion at NLO (33 fb) is also indicated (black dashed line) for comparison. One can see that the cross 
section for H^A^ is about 20-500 fb in the mass region 130-300 GeV, which is much larger than the hh 
production for most of the mass region and thus leads to great discovery potential. 

In general, the new scalars are nondegenerate for a nonzero A 5 . In the positive scenario where 
are the lightest, the cascade decays of and H^/A^ can strengthen the observability of ll67l l 6 ^ . 
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FIG. 9. Production cross sections for a pair of scalars at LHC14 versus Ma for a degenerate spectrum. The black 
dashed line is for the SM hh production. 

For the same reason, in the negative scenario where /A^ are the lightest, the charged scalars contribute 
instead to the production of H^/A^ through the cascade decays like —)• /A^W*. In this work, we 

study these contributions in the same way as was done for the positive scenario in Refs. llgTlIMl . 

We define the reference cross section Xq for the standard Drell-Yan process 

Xo = a{pp^ Z* ^ H^A^), (10) 

which is independent of the cascade decay parameters va. and AM. A detailed study on the bhT^T~ signal 
for this process with Ma = 130 GeV can be found in Ref. 1(1]. Besides the above direct production, neutral 
scalars can also be produced from cascade decays of charged scalars. These extra production channels 
include /A^, and followed by cascade decays of charged scalars. We 

consider first the associated H^H^/A^ production followed by cascade decays of 

pp^W* ^ pp^W* ^ A^H^W*, 

pp^W* ^ H^A^ H^A°W*, pp^W* ^ H^A^ A^A°W*, (11) 

resulting in three final states classified by a pair of neutral scalars: A^H^, H^H^, and A^A^. Noting that 
the last two originate only from cascade decays, any detection of such production channels would be a hint 
of charged scalars being involved. Using the fact that 


(j{pp ^ lU* ^ H^H^) ~ a{pp ^ FF* ^ H^A^) 
BR{H^ H^W*) ~ BR{H^ A°W*), 


( 12 ) 

( 13 ) 
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as well as the narrow width approximation, we calculate the production cross sections for these three final 
states: 

= 2[a{pp ^ ^ + a{pp ^ ^ R-H^)] x BR{H^ A'^W*), (14) 

: Yi = [a{pp ^ ^ H+H^) + a{pp ^ IV" ^ R-R^)] x BR{R^ R^W*), (15) 

A^A° : Zi = [a{pp ^ 14^+ ^ R+A^) + a{pp ^ ^ R-A°)] x BR{R^ A^W*). (16) 

The factor 2 in Xi accounts for the equal contribution from the process with R^ and AP interchanged. The 
relations Xi = 2Yi = 2 Zi actually hold true for all of the four production channels, since for a given 
channel the same branching ratios (such as for R^ R^/A^W*) are involved, 

Xi = 2Yi = 2Zi, {i = l,2.2,,A), (17) 

where Xi, Yi, and Zi refer to the cross sections for R^A^, R^R^, and A^A^ production with the subscript 
f = 1, 2, 3, 4 denoting the production channels R^R^/A^, R~^R~ , R^R^^, and R'^'^R , respectively. 

The relations imply that we may concentrate on the cross section of R^A9 production. 

Naively, one would expect the next important channel to be R^R~ since it only involves two cascade 
decays: 

X 2 = 2a{pp -i*lZ* R+R-) X BR(iT=^ ^ R^W*)BR{R^ A^W*). (18) 

But as already mentioned in Ref. ll68l . a smaller coupling and destructive interference between the 7 * and 
Z* exchange make the cross section of R^R~ production an order of magnitude smaller than that of R^AP 
even for a degenerate spectrum. Considering further suppression due to cascade decays, X 2 is not important 
for the enhancement of R^A9 production and can be safely neglected in the numerical analysis. 

The contribution from R^R^^ is more important despite the fact that it involves three cascade decays: 

X 3 = 2[a{pp W-* R+R—) + (t(pp ^ R-R++)] x (19) 

BR(i7=^=^ ^ R^W*)BR(R^ R°W*)BR(R^ A°W*). 

As shown in Fig.[^ a{pp —)■ W* —)• R^R^^) is the largest for a degenerate mass spectrum. When cascade 
decays are dominant, the phase-space suppression of heavy charged scalars will be important. So we expect 
that the R^A^ production receives considerable enhancement from R^R^^ when the mass splitting is 
small and cascade decays are dominant. 

Finally, the last mechanism is R^~^R , which involves four cascade decays: 

X 4 = 2a{pp -i*lZ* iT++iT“) X BR(iT=^=^ ^ R^W*f 
xBR(iT=^ ^ R°W*)BR{R^ A^W*) . 


( 20 ) 


14 




FIG. 10. Production cross sections for a pair of neutral scalars versus Ma at LHC14 and with AM = 5 GeV, 
UA = 0.001 GeV. Left: The red solid (dashed) line corresponds to Xq (X). Right: The red line corresponds to 
/A^A^ from cascade decays YjZ, and the green line to H^A^ from cascade decays {Xq)- The shaded regions 
are hlled by scanning over AM and va- 

This mechanism is also promising since the cross section of production is slightly larger than 

production for a degenerate mass spectrum. The phase-space suppression of X 4 is more severe than 
that of X 3 , because a pair of the heaviest are produced. 

Summing over all four of the above channels yields the contribution to the production from 

cascade decays, 


Xc = + X2 + X3 + X4, 


( 21 ) 


and the total production cross section of H^A^ is then X = Xq + Xq- Using Eq. (171, the total cross 
sections for the pair production /A^A^, Y = Yhi Yi, Z = Zi, are given by 


y = Z = -Xc. 
2 ^ 


( 22 ) 


Since the enhancement from cascade decays depends on a not severely suppressed phase space and a larger 
branching ratio of cascade decays, we choose to work with a relatively smaller mass splitting and triplet vev 
as shown in Eq. (j^. Eigure[^ displays the cross sections of the H^A^, H^H^, and A^A^ production as a 
function of Ma. As can be seen from the figure, the production of H^A^ can be enhanced by a factor of 
3, while the /A^A^ production at the maximal enhancement can reach the level of Xq. This could 

make the detection of neutral scalar pair productions very promising in the negative scenario. 
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IV. LHC SIGNATURES OF NEUTRAL SCALAR PRODUCTION 


In this section we investigate the signatures of neutral scalar production at the LHC. From previous 
studies on the SM hh production, we already know that the most promising signal is 6677 , and bhT~^T~ is 
next to it, while both semileptonic and dileptonic decays of FF’s in the hbW^W~ channel are challenging. 
In this work we analyze all three of the signals— 6677 , bbT~^T~, and bbW^W~ —)• bM^£~2iy {£ = e, ^ for 
collider identification)—as well as their backgrounds based on the benchmark model presented in Eq. Q. 

we discussed the Drell-Yan production of and the enhanced pair and associated pro- 


In Sec. 


Ill 


duction of neutral scalars H^/A^ due to cascade decays of charged scalars H^, We are now ready 

to incorporate the branching ratios of H^/A^ decays for a specific signal channel. For instance, the cross 
sections for the 6677 signal channel can be written as 

50 ( 6677 ) = Xo X [BR(ff° ^ 66 )BR(yl° ^ 77 ) + BR(ii'° ^ 77 )BR(yl° ^ 66 )] , (23) 

5 ( 6677 ) = ^ X [BR(ff° ^ 66)BR(2l° ^ 77 ) + BR(H'° ^ 77 )BR( 2 l° ^ 66 )] (24) 

+2yx BR{H° 66 )BR(iF° ^ 77 ) + 2Zx BR(^° ^ 66 )BR(A° ^ 77 ). 

Here So denotes the signal from the direct production pp ^ Z* ^ H^A^ alone, and 5 includes contribu¬ 
tions from cascade decays. 5(o)(66T'''r”) has a similar expression as 5 (o)( 6677 ), while is 

simpler since the decay mode A^ —)• W^W~ is absent. 

The theoretical cross sections for the 6677 , bbT~^T~, and bbl'^i~2v signal channels are plotted in 
The cross section So{blry^)/So{bbT~^T~) is larger than that of the SM hh production until 


Fig. 


11 


Ma = 159/161 GeV; taking into account cascade enhancement pushes the corresponding Ma further 
to 179/197 GeV. So{hbt^l~2u) is always larger than that of hh in the mass region 130-200 GeV, and 
interestingly, it keeps about the same value when Ma < 160 GeV. The signal from is comparable 

with So in these three channels only for Ma < 160 GeV, while in contrast the signal from A^A^ becomes 
dominant for the 6677 and bbT~^T~ channels when Ma > 160 GeV. Therefore, we have a chance to probe 


the A^A^ pair production in these two channels. Also shown in Fig. 11 is the enhancement factor S/So for 


the three signal channels at the benchmark point Q as a function of Ma, which will help us understand the 
simulation results. 


A. 6677 signal channel 

In our simulation, the parton-level signal and background events are generated with madgraphs Il69l . 
We perform parton shower and fast detector simulations with pythia ITTOII and delphes3 iTTll . Finally, 
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FIG. 11. Theoretical cross sections of 6677 , and bM~^£~2i/ signal channels at LHC 14. The red solid (dashed) 

line corresponds to the signal from Xq (X), the green (blue) solid line corresponds to the signal from Y {Z), and the 
purple dashed line shows the total cross section S for the signal. The SM hh cross section is shown for comparison. 
The lower right panel shows the enhancement factor S/Sq in the three signal channels. 


MADANALYSIS5 f/Hl is responsible for data analysis and plotting. We take a flat 6-tagging efficiency of 70%, 
and mistagging rates of 10% for c jets and 1% for light-flavor jets, respectively. Jet reconstruction is done 
using the anti-fc^ algorithm with a radius parameter of R = 0.5. We further assume a photon identification 
efficiency of 85% and a jef-faking-photon rate of 1.2 x 10“^ f/^ . 
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The main SM backgrounds to the signal are as follows: 

6677 : pp —>■ 6677 , (25) 

tth : pp —> tth —> M~iy 77 missed), (26) 

Zh : pp — Zh — 6677 . (27) 

Among them, 6677 and Zh are irreducible, while tih is reducible and can be reduced by vetoing the 
additional i’s with p^ > 20 GeV and \r]£\ < 2.4. In addition, there exist many reducible sources of fake 
6677 : 

pp —)• bbjj bb^^,pp —)• bbj^ ^ 6677,... 

pp —)• CC77 ^ 6677, pp —)• j j77 ^ 6677,. .., (28) 


where x ^ y stands for a final-state x misidentified as y. The remaining fake sources are subdominant 
and are thus not included in our simulation. The QCD corrections to the backgrounds are included by 
a multiplicative Tf-factor of 1.10 and 1.33 for the leading cross sections of tih and Zh at LHC14 ll74l . 
respectively. The cross section of the 6677 background has been normalized to include fake sources and 
does not take NLO corrections into account. 

The distributions of some kinematical variables before applying any cuts are shown in Fig. where 
we assume Ma = 130, 160, 190 GeV. In our analysis, we require that the final slates include exaclly one 
6 -jef pair and one 7 pair and salisfy Ihe following basic culs: 

p^^ > 30 GeV, \rjb,'y\ < 2.4, > 0.4, (29) 

where AR = \J(Acj))'^ + (At/)^ is Ihe particle separation, wifh A(p and Arj being fhe separation in fhe 
azimufhal angle and rapidify, respecfively. Here we employ a lighler px cul lhan is usually applied fo 
suppress fhe QCD-eleclroweak 6677 background. The 6 -jel pair and 7 pair are Ihen required lo fall in fhe 
following windows on fhe invarianf masses and fulfill fhe AR cul crileria: 


ARbb < 2.5, \Mbb — T7a| < 15 GeV, 
AR-y-y < 2.5, - MaI < 10 GeV. 


(30) 


As shown in Fig. 12 fhe ARbb^y'y dislribulions of fhe signal are clearly more compacl as Ihey are more 
likely coming from fhe same parlicles. Thus fhe AR cufs can effectively suppress fhe background. More 
specific culs are necessary for furlher analysis. A useful variable is fhe invarianf mass of fhe neufral scalar 
pair and fhe lolal Iransverse energy Et is also dislinclive. The peak of Muo^o increases wifh Ma, 

and similarly for Et- For simplicify, we adopf for fhe cufs a linear shiff belween , Et and Ma: 


MjjOy^o > 2Ma + 90 GeV, Et > 2Ma — 60 GeV. 


( 31 ) 
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FIG. 12. Distributions of for the signal 6677 and its backgrounds before 

applying any cuts at LHC14. 
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Ma = 130 GeV 

H°A°{So) 

6677 

tth 

Zh 

S/B 

S{S,B) 

Cross section at NLO 

8.01 X 10-1 

5.92 X 103 

1.18 

2.99 X 10-4 

1.39 X 10-4 

5.75 X 10-4 

Basic cuts 

1.22 X 10-1 

4.16 X lOi 

1.03 X 10-4 

3.41 X 10-2 

2.92 X 10-3 

1.03 

ReconsU'uct scalars from hs 

6.99 X 10-2 

7.07 

1.50 X 10-2 

9.61 X 10-4 

9.87 X 10-3 

1.44 

Reconstruct scalars from 7 s 

5.28 X 10-2 

1.03 X 10-4 

1.08 X 10-2 

7.32 X 10-4 

4.63 X 10-4 

8.01 

Cut on Mho AO 

4.21 X 10-2 

2.04 X 10-2 

4.69 X 10-3 

3.23 X 10-4 

1.65 

12.0 

Cut on Et 

3.31 X 10-2 

6.58 X 10-3 

4.68 X 10-3 

2.27 X 10-4 

2.88 

12.8 

Cascade enhanced 

1.51 X 10-1 

- 

- 

- 

13.1 

41.0 

Ma = 160 GeV 

H°A°{So) 

6677 

tih 

Zh 

S/B 

S{S,B) 

Cross section at NLO 

5.10 X 10-2 

5.92 X 103 

1.18 

2.99 X 10-4 

8.61 X 10-6 

3.63 X 10-2 

Basic cuts 

8.78 X 10-3 

4.16 X 104 

1.03 X 10-4 

3.41 X 10-2 

2.10 X 10-4 

7.44 X 10-2 

Reconstruct scalars from 6 s 

4.11 X 10-3 

5.06 

1.34 X 10-2 

2.36 X 10-4 

8.11 X 10-4 

9.99 X 10-2 

Reconstruct scalars from 7 s 

3.27 X 10-3 

3.42 X 10-2 

1.57 X 10-3 

0.00 

9.56 X 10-2 

9.53 X 10-4 

Cut on Mho AO 

2.57 X 10-3 

1.12 X 10-2 

1.18 X 10-3 

0.00 

2.30 X 10-4 

1.28 

Cut on Et 

1.73 X 10-3 

3.95 X 10-3 

1.03 X 10-3 

0.00 

4.37 X 10-4 

1.41 

Cascade enhanced 

1.10 X 10-2 

- 

- 

- 

2.77 

7.29 

Ma = 190 GeV 

H°A°{So) 

6677 

tth 

Zh 

S/B 

S{S,B) 

Cross section at NLO 

2.68 X 10-3 

5.92 X 103 

1.18 

2.99 X 10-4 

4.53 X 10-^ 

1.91 X 10-3 

Basic cuts 

5.33 X 10-4 

4.16 X 104 

1.03 X 10-4 

3.41 X 10-2 

1.28 X 10-3 

4.52 X 10-3 

ReconsU'uct scalars from 6 s 

2.27 X 10-4 

3.61 

1.05 X 10-2 

1.24 X 10-4 

6.27 X 10-3 

6.53 X 10-3 

Reconstruct scalars from 7 s 

1.81 X 10-4 

2.47 X 10-2 

3.93 X 10-6 

0.00 

7.34 X 10-3 

6.30 X 10-2 

Cut on Mho AO 

1.55 X 10-4 

9.87 X 10-3 

3.93 X 10-6 

0.00 

1.57 X 10-2 

8.52 X 10-2 

Cut on Et 

8.35 X 10-3 

1.48 X 10-3 

3.93 X 10-6 

0.00 

5.63 X 10-2 

1.18 X 10-4 

Cascade enhanced 

1.50 X 10-3 

- 

- 

- 

1.01 

1.87 


TABLE I. Evolution of signal and background cross sections (in fb) at LHC 14 for the blr/j signal channel upon 
imposing the cuts one by one. Eor the cascade-enhanced signal only the cross section passing all the cuts is shown. 
The last two columns assume an integrated luminosity of 3000 fb^^. 

For instance, we apply > 350 GeV, Et > 200 GeV at the benchmark point Ma = 130 GeV. 

To estimate the observability quantitatively, we adopt the following significance measurement: 

-S-C^ (32) 

which is more suitable than the usual definition of S/^/B or S/\/S + B for Monte Carlo analysis fTSll . 
Here S and B are the signal and background cross sections, and C is the integrated luminosity. The survival 
cross sections of the signal from the Drell-Yan process and of the backgrounds upon imposing cuts step by 


5(5,H) = ^/2((S-£ + H.£)log(l + | 
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step are summarized in Tablej^at the benchmark point Q for Ma = 130, 160, 190 GeV respectively. For 
the cascade-enhanced signal, only the cross section passing all the cuts is shown. The last two columns in 
the table show the signal-to-background ratio S/B and the statistical significance S{S, B). 

For Ma = 130 GeV, the 6677 channel is very promising. Without (with) cascade enhancement, the final 
significance can reach 12.8 (41) for LHC 14@3000, corresponding fo 99 (453) evenfs. For Ma = 160 GeV, 
fhe channel becomes challenging since fhe cross secfion has decreased by a factor of 15.7 compared wifh 
fhe case of Ma = 130 GeV. Buf fhe cufs we applied are efficienl fo suppress fhe SM background, and wifh 
cascade enhancemenf fhe significance could sfill reach 7.29 for 3000 fb~^, corresponding to 33 evenfs in 
fhe mosf optimistic case. For Ma = 190 GeV, if looks hopeless even wifh maximal cascade enhancemenf 
in our benchmark model: fo achieve 10 signal evenfs, an integrated luminosify of af leasf 6670 fh“^ is 
required, which is beyond fhe reach of fhe fufure LHC. 


B. bbr^T signal channel 


For fhis signal channel, an imporfanf pari of fhe analysis depends on fhe abilify fo reconslrucf fhe b pair 
and fhe r pair. Here we consider fhe hadronic decays of fhe r leplon and assume a r-lagging efficiency of 
70% wifh a negligible fake rale. 

The main SM backgrounds are as follows: 


bbT~^T~ : pp —)■ bbZ/j*/h —)• bbT^T~, 
bbW~^W~ : pp —)• bbW^W~ —)• bbT^VrT~i'T, 
Zh : pp —)■ Zh —)■ bbT~^T~. 


(33) 

(34) 

(35) 


The irreducible QCD-elecfroweak background comes from bbT^T~, where fhe r pair originales from fhe 
decays of Z/7*//i. Since fhe hadronic decays of r always confain neufrinos, we also include fhe SM back¬ 
ground bbW~^W~, which confribules to fhe hbT~^VrT~i'T final slate. The hbW~^W~ background mainly 
originales from tt produclion wifh subsequenl decays t —)• bW and W —)■ TVr- Moreover, fhe associated 
Zh production gels involved Ihrough fhe subsequenl decays h ^ bb and Z —)■ or vice versa. The 

QCD corrections to fhe backgrounds are included by a multiplicative iT-faclor of 1.21, 1.35, and 1.33 to fhe 
leading cross section of ftfer+r” fT^ . tt IfTTl . and Zh llTdll af LHC14. 


The kinemalical dislribulions similar to fhe 6677 channel are shown in Fig. 13 As one can see from fhe 
figure, fhe r jefs are less energetic lhan fhe 6 jels (similar fo Ihose in fhe 6677 signal channel) due fo missing 
neufrinos in fhe final sfafe. We firsl employ fhe following selection cufs fo pick up signals wifh exacfly one 
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FIG. 13. Distributions of ARtt, Et, andfor the signal bbr^r and its backgrounds before 

applying any cuts at LHC14. 


b pair and one r pair: 

p'rf > 30 GeV, \r]b^r\ < 2.4, ARbb,bT,TT > 0.4, (36) 

and no cut on is adopted. After the selection, the r and b pairs are required to fulfill the cuts on the 
invariant masses and separations: 


ARrr < 2.5, Ma - 40 GeV < Mrr < Ma 
AR bb < 2.5, \Mbb — Ma\ < 15 GeV. 


( 37 ) 
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Ma = 130 GeV 

H^A^iSo) 

bbr'^T 

bbW+W- 

Zh 

S/B 

SiS,B) 

Cross section at NLO 

4.31 X 10^ 

3.10 X 104 

7.92 X 103 

2.21 X IQi 

1.11 X 10-3 

12.0 

Basic cuts 

7.75 X 10-1 

4.49 X IQi 

8.97 X lOi 

2.91 X 10-1 

5.74 X 10-3 

3.65 

Reconstruct scalars from rs 

5.14 X 10-1 

1.19 X IQi 

3.57 X lOi 

1.06 X 10-1 

1.08 X 10-2 

4.06 

Reconstruct scalars from 6s 

2.14 X 10-1 

4.34 

9.44 X 10-1 

2.28 X 10-2 

4.04 X 10-2 

5.06 

Cut on MffOAo 

1.29 X 10-1 

1.96 X 10-1 

1.51 X 10-1 

8.10 X 10-3 

3.64 X 10-1 

11.3 

Cut on Et 

1.03 X 10-1 

9.87 X 10-2 

7.35 X 10-2 

5.89 X 10-3 

5.81 X 10-1 

12.4 

Cascade enhanced 

5.27 X 10-1 

- 

- 

- 

2.96 

51.6 

Ma = 160 GeV 

H^A°{So) 

bbT'^T~ 

bbW+W- 

Zh 

S/B 

S{S,B) 

Cross section at NLO 

3.08 

3.10 X 104 

7.92 X 103 

2.21 X IQi 

7.92 X 10-3 

8.55 X 10-1 

Basic cuts 

6.81 X 10-2 

4.49 X IQi 

8.97 X lOi 

2.91 X 10-1 

5.05 X 10-4 

3.21 X 10-1 

Reconstruct scalars from rs 

3.14 X 10-2 

1.52 X IQi 

2.46 X 102 

3.17 X 10-2 

1.20 X 10-3 

3.36 X 10-1 

Reconstruct scalars from 6s 

1.2 X 10-2 

2.47 

1.06 X 10-1 

0.00 

4.80 X 10-3 

4.10 X 10-1 

Cut on MffOAo 

6.99 X 10-3 

1.22 X 10-1 

2.06 X 10-2 

0.00 

4.89 X 10-2 

1.00 

Cut on Et 

5.04 X 10-3 

4.72 X 10-2 

5.88 X 10-3 

0.00 

9.48 X 10-2 

1.18 

Cascade enhanced 

5.11 X 10-2 

- 

- 

- 

9.63 X 10-1 

10.7 

Ma = 190 GeV 

H°A°{So) 

bbT'^T~ 

bbW+W- 

Zh 

S/B 

S{S,B) 

Cross section at NLO 

2.47 X 10-1 

3.10 X 104 

7.92 X 103 

2.21 X lOi 

6.34 X 10-3 

6.86 X 10-2 

Basic cuts 

6.54 X 10-3 

4.49 X IQi 

8.97 X lOi 

2.91 X 10-1 

4.86 X 10-3 

3.09 X 10-2 

Reconstruct scalars from rs 

2.32 X 10-3 

4.60 X 10-1 

1.47 X lOi 

5.89 X 10-3 

1.53 X 10-4 

3.26 X 10-2 

Reconstruct scalars from 6s 

7.66 X 10-4 

2.06 X 10-2 

1.21 

0.00 

6.25 X 10-4 

3.78 X 10-2 

Cut on MffOAo 

6.34 X 10-4 

1.47 X 10-3 

9.03 X 10-2 

0.00 

6.91 X 10-3 

1.14 X 10-1 

Cut on Et 

3.87 X 10-4 

0.00 

2.64 X 10-2 

0.00 

1.47 X 10-2 

1.30 X 10-1 

Cascade enhanced 

6.85 X 10-3 

- 

- 

- 

2.60 X 10-1 

2.22 


TABLE n. Similar to Tableji] but for the bbr^r signal channel. 

The different mass shift between M-tt and Mbb is owing to the missing neutrinos in r decays resulting in a 
wider distribution of Mrr- For the reconstructed neutral scalars, we further adopt similar cuts on M^oao 
and Et as in the 6677 channel: 

Mroao > 2 M/^ + 70 GeV, Ex > 2 M/^ — 80 GeV. (38) 

Both Mjjoao and Et cuts are reduced by 20 GeV compared with the 6677 channel, which again results 
from neutrinos in the final state. The corresponding results are summarized in Table |n| 

The fthr+r” is also promising for Ma = 130 GeV even without enhancement from cascade decays. 
The final significance is 12.4 and the corresponding number of signal events is 309 for LHC 14@3000. 
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Including the cascade enhancement, the significance is improved to 51.6, which is even better than the 6677 
signal. For Ma = 160 GeV, the biggest challenge is also the small production cross section of the signal. 
But in the most optimistic case, the cascade decays can increase the signal by a factor of 10.1, making 
this channel feasible. Finally, neutral scalars as heavy as 190 GeV are difficult to detect at LHC 14 in this 
channel. 


C. bbW^W signal channel 


It is difficult to search for the SM Higgs pair production in this channel due to missing energy brought 
about by neutrinos in leptonic decays of the W boson, which makes one of the two Higgs bosons not fully 
reconstructible |[3^l57l . The situation is ameliorated in our scenario because, the production rate of 
can be an order of magnitude larger than that of hh and the Ai-W decay branching ratio of can also 
be larger than h in the vast parameter space. This considerably improves the signal events and partially 
compensates the deficiency of the detection capability. 

With both VF’s decaying leptonically, the final state appears as bbi~^£~^T- The dominant SM back¬ 
grounds are as follows: 

ti :pp^ti^ bW+bW- bU+r^r- (39) 


As before, the QCD correction is included by a multiplicative iT-factor of 1.35 for the ti production fTTl . 
We pick up the events that include exactly one 6 -jet pair and one opposite-sign lepton pair and filter them 
with the basic cuts: 


> 30 GeV, p^ > 20 GeV, < 2.4, 
^Rbb,MM > 0.4, Xr > 20 GeV. 

The separation and invariant mass of the 6 -jet pair are required to fulfill 

< 2.5, iMfefe - MaI < 15 GeV. 

For the lepton pair, we reconstruct the transverse cluster mass Mq\ 


(40) 


= 


\l {\IPT,ie + + [pt,u +Xr) 


(41) 


(42) 


The distributions of M^, ARu, andare shown in Fig. 14 The peak of is always lower than Ma by 
about 30-40 GeV, and the lepton separation AR^^ in the signal is much smaller than in the ti background. 
Accordingly, we set a wide window on while tightening up the cuts on AR^ and 


Ma ~ 80 GeV < M^ < Ma, AR^ < 1 . 2 , > 0.9Ma. 


(43) 
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FIG. 14. Distributions of ARu, Mq, ^t, Mq, and Et for the signal bbi'^i 0 t and its backgrounds before 
applying any cuts at LHC14. 

We find that is least efficient around Ma ~ 190 GeV, where the peak of for the ti background 
is around 150 GeV. The very tight cuts on AR^i and are sufficient to suppress the background by 1 or 
2 orders of magnitude, while keeping the number of signal events as large as possible. We further combine 
the 6 -jet pair and the lepton pair into a cluster and construct the transverse cluster mass: 

Me = y ( PT,bbu + +Xr) - {pT,bbee , (44) 

which is an analog of Mj^o^o in the previous subsection. The distribution of Me is displayed in Fig. [T^ 
which is very similar to that of in the 6677 channel. Although it looks from the Me distributions 
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Ma = 130 GeV 


tt 

S/B 

S{S,B) 

Cross section at NLO 

3.91 

2.38 X lOi 

1.69 X 10-4 

1.41 

Basic cuts 

1.51 

4.04 X 103 

3.74 X 10-4 

1.30 

Reconstruct scalars from 6s 

3.29 X 10-1 

3.35 X 102 

9.82 X 10-4 

0.984 

Cut on Mq 

3.21 X 10-1 

2.14 X 102 

1.50 X 10-3 

1.20 

Cut on ARii 

2.64 X 10-1 

9.26 X lOi 

2.85 X 10-3 

1.50 

Cut on 

8.45 X 10-2 

1.48 X lOi 

5.71 X 10-3 

1.20 

Cut on Me 

3.30 X 10-2 

1.69 X 10-1 

1.95 X 10-1 

4.26 

Cut on Et 

3.19 X 10-2 

1.47 X 10-1 

2.17 X 10-1 

4.41 

Cascade enhanced 

1.40 X 10-1 

- 

9.53 X 10-1 

17.7 

Ma = 160 GeV 

H^A^iSo) 

ti 

S/B 

S{S,B) 

Cross section at NLO 

4.95 

2.38 X lOi 

2.13 X 10-4 

1.78 

Basic cuts 

2.13 

4.04 X 103 

5.27 X 10-4 

1.84 

Reconstruct scalars from 6s 

4.25 X 10-1 

2.68 X 102 

1.59 X 10-3 

1.42 

Cut on 

3.97 X 10-1 

1.89 X 102 

2.10 X 10-3 

1.58 

Cut on ARii 

3.21 X 10-1 

7.04 X lOi 

4.56 X 10-3 

2.09 

Cut on 

9.47 X 10-2 

4.29 

2.21 X 10-2 

2.50 

Cut on Me 

3.28 X 10-2 

4.74 X 10-2 

6.92 X 10-1 

7.50 

Cut on Et 

3.02 X 10-2 

3.62 X 10-2 

8.34 X 10-1 

7.78 

Cascade enhanced 

1.01 X 10-1 

- 

3.24 

23.2 

Ma = 190 GeV 

H^A^iSo) 

ti 

S/B 

S{S,B) 

Cross section at NLO 

1.19 

2.38 X lOi 

5.00 X 10-3 

0.424 

Basic cuts 

6.44 X 10-1 

4.04 X 103 

1.59 X 10-4 

0.554 

Reconstruct scalars from 6s 

1.36 X 10-1 

2.26 X 102 

6.02 X 10-4 

0.495 

Cut on Mq 

1.27 X 10-1 

1.79 X 102 

7.09 X 10-4 

0.520 

Cut on ARii 

9.70 X 10-2 

6.05 X lOi 

1.60 X 10-3 

0.683 

Cut on ^T 

2.57 X 10-2 

1.62 

1.59 X 10-2 

1.10 

Cut on Me 

8.85 X 10-3 

1.89 X 10-2 

4.68 X 10-1 

3.29 

Cut on Et 

8.37 X 10-3 

1.40 X 10-2 

5.98 X 10-1 

3.56 

Cascade enhanced 

2.69 X 10-2 

- 

1.92 

10.1 


TABLE III. Similar to Table|ij but for the bb£~^£ signal channel. 


(before any cuts are made) that the tt background has a large overlap with the signal, the cuts on M^, 
ARgi, andactually modify them remarkably, so that a further cut on Me could improve the significance 
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MaCGcV) MA(GeV) 


FIG. 15. Left; Significance S{S,B) of the 6677 channel versus Ma reachable at LHC14@300 (red region) and 
LHC14@3000 (green). Right; Required luminosity to reach a Scr (red region) and her (green) significance in the 
6677 channel versus Ma at LHC 14. The solid line corresponds to the signal from Xq alone, and the dashed line 
corresponds to the total signal including cascade enhancement. 


efficiently. We apply a cut on Me as we did with as well as one on Et'- 


Me > 2Ma + 90 GeV, Et > 2Ma — 60 GeV. 


(45) 


The results following the cutflow are summarized in Table III For Ma = 130 GeV, the final significance 


is 4.41 (17.7) without (with) cascade enhancement. With cascade enhancement this should be enough to 
discover the neutral scalars. The signal channel is more promising for Ma = 160 GeV due to a slightly 
larger cross section and higher cut efficiencies. The final significance is 7.78 (23.2), which is also better 
than the 6677 and hbT~^T~ channels with the same mass. Finally, for Ma = 190 GeV, the significance 
becomes 3.56 (10.1). Therefore, for our benchmark model, the only promising signal for such heavy neutral 
scalars (~ 190 GeV) comes from the bhW~^W~ channel. 


D. Observability 


Based on our elaborate analysis of signal channels in Secs. IV A IVC[ we examine the observability 
of the neutral scalars Mq, Aq in the mass region 130 ~ 200 GeV by adopting essentially the same cuts 


as before. In the left panel of Figs. 15 16 and 17 we present the significance S{S, B) as a function of 


Ma in the three signal channels 6677, bbT~^T , and bbt^l that is reachable for LHC14@300 and 
LHC14@3000, respectively. The required luminosity to achieve a 3a and 5a significance is displayed in 
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FIG. 16. Same as Fig. 15 but for the bbr^r channel. 


the right panel of the figures. As was done in our previous analysis, the effect of cascade enhancement is 
included by a factor S/S'o in the final resulfs. 


As shown in Figs. 15 and 16 bofh fhe 6677 and fthr+r” channel are fypically sensifive fo fhe low-mass 
region (Ma < 160 GeV). In fhe absence of cascade enhancemenf, fhe 3a significance would never be 
reached for M/\ > 138 (142) GeV in fhe 6677 {bbT^T~) channel for LHC14@300. However, a cascade 
enhancement of 5/S'o ~ 4 — 6 (as can be seen from Fig. [IT]) in this mass region can greatly improve the 
observability, pushing the 3a mass limit up to 157 (162) GeV in the 6677 (bbT~^T~) channel. Moreover, 
with cascade enhancement, one has a good chance to reach a 5a significance if Ma < 153 (155) GeV. In 
ofher words, fhe cascade enhancemenf significanfly reduces fhe required luminosify. For insfance, fo achieve 
a 3cr and 5a significance in fhe 6677 (bbT~^T~) channel wifh Ma = 130 GeV, fhe required luminosify is 
as low as 16 (10) fb~^ and 42 (27) fb“^ af LHC14, respecfively. The channel is more promising, 

fhanks fo a relafively larger production rafe. 

Af fhe fufure LHC 14 wifh 3000 fb“^ dafa, fhe heavier mass region can also be probed. Wifh a maximal 
cascade enhancemenf, fhe 3 cj and 5a mass reach is pushed fo 177 and 164 GeV, respecfively, in fhe 6677 
channel, which should be compared fo 156 and 151 GeV in fhe absence of enhancemenf. For fhe bbT~^T~ 
channel, fhe enhancemenf facfor S/Sq can reach abouf 18 above fhe VF-pair fhreshold, upshiffing fhe 3a 
and 5cj mass reach fo 189 and 177 GeV, respecfively, from 154 and 150 GeV wifhouf fhe enhancemenf. 

The bM^£~^T channel shown in Fig. 17 is more special, compared wifh 6677 and bbT^T~. If is 
relafively more sensifive fo a higher mass befween 150-180 GeV, where fhe decay mode W~^W~ 

dominafes, while ifs observabilify deferiorafes for Ma < 150 GeV due fo phase-space suppression in fhe 
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FIG. 17. Same as Fig. 15 but for the bM^i channel. 


decay. The cascade enhancement S/Sq at our benchmark point (j^ is typically 3-4 in the mass region 130- 
200 GeV, and decreases as Ma increases. For LHC14@300, the 3a and 5 ct mass reach is, respectively, 
190 and 181 GeV with maximal cascade enhancement. These limits would just increase by 2-3 GeV for 
LHC14@3000 if there were no cascade enhancement, while with cascade enhancement the 5a limit, for 
instance, is pushed up to 200 GeV. Finally, a 3fT or 5cr reach in the mass region 150-180 GeV requires an 
integrated luminosity of 50 fb~^ (450 fb“^) or 150 fb“^ (1300 fb“^) with (without) cascade enhancement. 


V. DISCUSSIONS AND CONCLUSIONS 

In this paper, we have systematically investigated the LHC phenomenology of neutral scalar pair pro¬ 
duction in the negative scenario of the type II seesaw model. To achieve this goal, we first examined the 
decay properties of the neutral scalars Hq/Aq and found that the scalar self-couplings A* have a great impact 
on the branching ratios of /A^. The coupling A 4 is important for tree-level decays of and AP, while 
one-loop-induced decays of further depend on A 2 and A 3 . We found that the decay —)• W~^W~ 
could dominate for 2Mw < Mjjo < 2Mh with A 4 < 0, while it can be neglected once Mj^o is above the 
light scalar pair threshold 2M/j. Moreover, the branching ratios of the decays —)■ 77 , Z 7 can cross 3 
orders of magnitude when varying the couplings A*, and there exist zero points for the H^ZZ, H^hh, and 
APZh couplings. 

The cross section of the Drell-Yan process pp ^ Z* ^ H^A^ for Ma < 200 GeV is much larger than 
that of the SM Higgs pair production driven by gluon fusion. In this paper, we studied the contributions 
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to production from cascade decays of the charged scalars and There are actually three 

different states for the neutral scalar pair: H^A^, H^H^, and A^A^. Here, and A^A^ can only 

arise from cascade decays of charged scalars, and their production rates always stay the same to a good 
approximation. Further, for a fixed value of Ma, cascade enhancement is determined by the variables 
VA and AM. By tuning these two variables, the associated production rate of H^A^ can be maximally 
enhanced by about a factor of 3, while those of the and A^A^ pair production can reach the value of 

H^j \0 production through the pure Drell-Yan process. 

We implemented detailed collider simulations of the associated H^A^ production for three typical signal 
channels bbT^T~, and bbW~^W~ with both W’s decaying leptonically). The enhancement from 

cascade decays of charged scalars is quantified by a mulfiplicafive facfor S/So. Due mainly fo a larger 
production rafe, all fhree channels are more promising fhan fhe SM Higgs pair case. If fhere were no cascade 
enhancemenf, fhe 5 cj mass reach of fhe 6677 , and bM^i~^T channels would be, respectively, 151, 

150, and 180 GeV for LHC14@3000. The cascade enhancemenf pushes fhese limifs up fo 164, 177, and 
200 GeV. The 6677 and bbT~^T~ channels are more promising in fhe mass region below abouf 150 GeV, 
and fhe required luminosifies for 5 cj significance are 42 fb~^ and 27 fb“^, respectively, af our benchmark 
poinf. Compared wifh fhese fwo channels, fhe bM~^£~^T channel is more advanfageous in fhe relafively 
higher mass region 150-200 GeV, and fhe required luminosify for 5 cj significance is abouf 150 fb“^ wifh 
maximal cascade enhancemenf. Needless fo say, for fhe purpose of a full invesfigafion on fhe impacf of 
heavy neufral scalars on fhe SM Higgs pair producfion, more sophisficafed simulafions are necessary. We 
hope fhaf fhis work may shed some lighf on furfher sfudies in bofh fhe phenomenological and experimenfal 
communities. 
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